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We investigated the nature of the magnetic phase transition in Gei-aj-ySn^MnyTe mixed crystals 
with chemical composition changing in the range of 0.083 < x <0. 142 and 0.012 < y < 0.119. The 
DC magnetization measurements performed in magnetic fields up to 90 kOe and temperature range 
2 -^ 200 K showed that the magnetic ordering at temperatures below T = 50 K exhibits features 
characteristic for both spin-glass and ferromagnetic phases. The modified Sherrington - Southern 
model was applied to explain the observed transition temperatures. The calculations showed that 
the spin-glass state is preferred in the range of the experimental carrier concentrations and Mn 
contents. The value of the Mn hole exchange integral was estimated to be Jpd = 0.45±0.05 eV. 
The experimental magnetization vs temperature curves were reproduced satisfactorily using the 
non-interacting spin-wave theory with the exchange constant Jpd values consistent with those cal- 
culated using modified Sherrington - Southern model. The magnetization vs magnetic field curves 
showed nonsaturating behavior at magnetic fields i? < 90 kOe indicating the presence of strong 
magnetic frustration in the system. The experimental results were reproduced theoretically with 
good accuracy using the molecular field approximation-based model of a disordered ferromagnet 
with long-range RKKY interaction. 

PACS numbers: 72.80.Ga, 75.40.Cx, 75.40.Mg, 75.50.Pp 



I. INTRODUCTION 

The IV- VI group of ferromagnetic compounds e.g. 
Gei-2;Mn2;Te possesses unique possibility to tune the 
magnetic and electrical properties of the crystals in- 
dependently [1[. The magnetism of bulk (Ge,Mn)Te 
crystals was for the first time studied over 40 years 
ago by Rodot et. al 0. The ferromagnetic ordering 
observed in the Gci-^^Mn^^Te with the Curie temper- 
atures Tc as high as 167 K [3J for the crystals with 
the composition a; = 0.5 was attributed to the indirect 
long range Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction. Recently, molecular beam epitaxy (MBE) 
grown Gei-a;Mna;Te thin layers are the subject of great 
scientific interest^"—. The progress in the methods of 
growth of this compound allowed the increase of the 
observed Curie temperatures to values as high as 200 K 
[y] and gives hope to obtain the room temperature 
ferromagnetism. 

In this paper we extend our previous 
studies of dynamic mag netic properties of 
Gei-x-ySn^MnyTe crystals [13, UM- The results of 
the measurements of the AC susceptibilities showed 
the occurrence of the magnetic transitions into the 
spin- glasslike state at temperatures lower than 50 K. In 



this paper we present the measurements of the static 
magnetic properties i.e. the magnetization measure- 
ments in the presence of a static magnetic field. 

The observed magnetic properties are discussed in 
theoretical context. The Sherrington-Southern model, 
capturing the essentials of spin-glass physics, allows 
us to predict the relevant transition temperatures and 
to estimate the value of the exchange integral Jpd 
in the studied crystals, assuming the RKKY interac- 
tion between magnetic impurities. The temperature 
dependence of magnetization will be modeled using 
the noninteracting spin wave theory, including the 
same RKKY interaction, which supports the estimates 
of Jpd- To describe the behavior of the high- field 
magnetization curves, we construct a molecular field 
approximation-based model of a disordered ferromagnet 
with long-range RKKY interaction, putting special 
emphasis on the possible magnetic inhomogenity in 
the system, caused by presence of antifcrromagnetic 
superexchange interaction between nearest-neighbor 
magnetic ions. This model is used to explain the 
non-saturating behavior of magnetization. 
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II. SAMPLE PREPARATION AND BASIC 
CHARACTERIZATION 

For purposes of this research bulk Gei-a;-j,Sna;Mnj,Te 
crystals were prepared using modified Bridgman method. 
The growth method was modified according to the 
ideas proposed by Aust and Chalmers for the alumina 
crystals^^. The radial temperature gradient was present 
in the growth furnace, creating 15 deg slope of the crys- 
tallization front. It allowed the decrease of the crystal 
blocks in the ingot from a few down to a single one. 

The chemical composition of the alloy was determined 
using the x-ray fluorescence method (with typical relative 
uncertainty of about 10%). The as grown ingots were 
cut into 1 mm thick slices perpendicular to the growth 
direction. The measured chemical content changed con- 
tinuously along the growth direction. For the present in- 
vestigations we used the very same samples as presented 
in Ref. pjj i.e. the studied crystals had chemical com- 
position changing in a range of values 0.083 <x< 0.142 
and 0.012 <y< 0.119. The crystal slices selected for the 
present investigations had small inhomogeneity of the 
chemical composition of the order of the uncertainty of 
x-ray fluorescence method i.e. less than 10% of the re- 
sultant molar fraction x 

The crystallographic quality of Gei-a;-ySna;Mnj,Te crys- 
tals was studied by means of the x-ray diffraction (XRD). 
Analysis of the XRD measurement results showed that 
the studied crystals were single phased. They crystallized 
in the [111] direction distorted NaCl i.e. the rhombohe- 
dral structure. The lattice parameters obtained for the 
studied crystals had values similar to those reported in 
the literature for germanium telluride bulk crystals e.g. 
the lattice parameter a was around 5.98 A and the angle 
of distortion a was around 88.3 [deg] |13j . 

Basic electrical properties of Gei-ai-ySna^Mn^Te crys- 
tals were also characterized. The measurements (using 
the standard DC six contact technique) of the resistivity 
parallel to the current direction p^x and the Hall effect 
(using constant magnetic field B — IA kOe) for each sam- 
ple at room temperature were performed. The samples 
had parallelepiped shape with typical diameters around 
1x1x8 mm. The results showed that the studied crystals 
were highly degenerated p-type semiconductors with high 
hole concentration n > 10^^ cm~^ and low carrier mobil- 
ities /x< 100 cm^/(Vs). For more detailed discussion of 
electrical properties of Gei-^^-ySnajMuyTe see Refs..lOiand 

[m. 



III. RESULTS AND DISCUSSION 

The DC magnetometry was used in the present studies 
of the Gei-2:-j,Sna;Mnj,Te crystals. The very same sample 
pieces as in the case of transport characterization with 
electrical contacts removed were studied by means of DC 
magnetometry. The magnetization measurements were 
performed using two magnetometers. The temperature 



dependent magnetization in a small magnetic fields was 
collected using a Quantum Design SQUID MPMS XL- 
5 magnetometer. The isothermal magnetization curves 
and hysteresis loops were measured using the extrac- 
tion method employed by the LakeShore 7229 DC mag- 
netometer. This method has the lower accuracy than 
SQUID magnetometer, but allowed us to measure at 
magnetic fields as high as B = 90 kOe. 



III.l. Low field magnetization 

The temperature dependencies of the zero-field-cooled 
(ZFC) and field-cooled (FC) magnetization were mea- 
sured in the temperature range between 4.3 and 100 K 
using the constant magnetic field B — bQ Oe. Typical re- 
sults of the temperature dependencies of the magnetiza- 
tion for selected Gei-a;-ySna;Mnj,Te crystals with constant 
Sn and Mn contents varied between 0.039 and 0.115 are 
presented in Figure [TJ The differences between the FC 
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FIG. 1. Experimental (points) and theoretical (lines) mag- 
netization as a function of temperature for samples cooled in 
the absence (closed symbols) and the presence (open sym- 
bols) of external magnetic field _B = 50 Oe for two selected 
Gei-aj-j/SnajMuyTe samples with different chemical composi- 
tion (labels). 

and ZFC curves for the studied crystals were observed. 
Bifurcations between ZFC and FC magnetization values 
at low temperatures are typical features of the appear- 
ance of a spin-glass phase. 

The paramagnet-spin-glass transition temperatures 
TsG were determined from the temperature dependent 
magnetization curves at the bifurcation points. The val- 
ues obtained in this way are in good agreement with our 
previously determined Tsg values from AC magnetom- 
etry results (see Ref. [ill) . Both the dynamic and static 



magnetic properties of Gei-a;-j,Sna;Mnj,Te crystals show 
features similar to those observed in canonical spin-glass 
systems e.g. CuMn or AuFe diluted metallic alloys^^. 

The presented experimental FC M{T) curves for weak 
field were theoretically reproduced with the use of the 
non-interacting spin wave theory. The magnon disper- 
sion relation has been assumed in the isotropic form 

I 



e (/c) = Dfc^, with spin- wave stiffness coefficient D = 
^yS ^j J [Rij) R^p where J {Rij) is the RKKY exchange 
integral, given by a formula (H]), calculated for a pair of 



ions in the lattice with a mutual distance of Ri 



In 



the present considerations it is assumed that the indirect 
RKKY interaction with the exchange constant Jrkky 
has the following forroi^"— 



J RKKY {Rij) — N\ 



327r3n2 



sm{2kpRij) — 2kpRij cos{2kpRij) 
(2fcF%)4 



exp( — j^), 



(1) 



where tti* is the effective mass of the carriers, Jpd is the 
exchange integral between conducting holes and Mn ions, 
oo is the lattice parameter, kp = {^ir'^n/hNv) is the 
Fermi wave vector, Ny is the number of the valleys of 
the valence band, h is the Planck constant divided per 
27r, and Rij is the distance between magnetic ions. For 
nearest-neighbor ions, the additional antiferromagnetic 
superexchange coupling, parameterized by Jaf, has been 
included. This results from averaging the magnon en- 
ergies over possible orientations of wavevector k for its 
given length k (similar to the calculation in Ref. |j). 

We emphasize that using such a direction-averaged dis- 
persion relation for our system is motivated by the pres- 
ence of spatial disorder in a diluted magnetic system, 
which makes the influence of the exact lattice structure 
(involved in calculating the Fourier transform) less pro- 
nounced. Moreover, we limit our considerations to a 
quadratic part of the dispersion relation instead of us- 
ing the full expression valid for the whole Brillouin zone. 
This reflects the fact that in a strongly diluted system the 
typical distances between the magnetic ions are notice- 
ably larger than the lattice constant itself, so that only 
low-energy spin waves (of wavelength exceeding the men- 
tioned distance) constitute well-defined excitations. The 
mentioned assumptions lead tho the temperature depen- 
dence of magnetization following exactly the Bloch's law 
in the form 



Af(r) = Af(0) 1-{T/Tsw) 



^3/2 



(2) 



The characteristic temperature Tsw can be expressed as 



Tsw 



27rDa2 



C(3/2)j ■ 



(3) 



where a is the lattice constant and C, (3/2) ~ 2.612 is 
the appropriate Riemann zeta function. The performed 
calculations allowed us to satisfactorily reproduce the 
experimentally observed M{T) curves (see Fig. [IJ for 
the Jpd exchange coupling constant values consistent 
with those obtained from the further analysis of Curie- 
Weiss and spin-glass transition temperatures within the 
Sherrington-Southern model. It is worth mentioning that 
such a robust T^^^-like behavior of M{T) curves has also 



been noticed by 18] and |T9| . 

The proposed theoretical approach was not able to take 
fully into account the magnetic frustration in the studied 
system, which is reflected in the shape (long tail) of the 
M{T) curve at T^Tsg- 



III. 2. Estimation of the exchange integral Jj 



pd 



The spin-glass state in the diluted magnetic material 
is well described by the model proposed by Sherrington 
and Southern's. The Sherrington-Southern (SS) model 
was found to be the most appropriate approach for the 
description of systems with a large number of neighbor- 
ing paramagnetic ions. Moreover, the modified SS model 
was applied to compounds belonging to the group of IV- 
VI semimagnetic semiconductors i.e. Sni-a;Mna;Te and 
Phi-x-ySuxMuyTe [2l|. In the present work we adapted 
the modified SS model 21] to quantify the Jpd exchange 
interaction in the studied Gei-j;-ySn3;MnyTe mixed crys- 
tals. Sherrington and Southern proposed the Gaus- 
sian distribution of the interaction strength in Heisen- 
berg Hamiltonian, with mean value Jq and width AJ. 
The model allows the calculations of both paramagnetic 
Curie- Weiss O and paramagnet-spin-glass transition Tsg 
temperatures given by Eqs.[3]and[S] 



2S{S+l)x 
^- 3kn ^°' 

rp 2AJ|- 2^C , 1\2 , O^O : 1^/0^1/2 



fSG 



3fcf 



■[S^iS + if + 8(8+ l)/2] 



(4) 
(5) 



where ks is the Boltzmann constant. In the above 
equations the RKKY indirect interactions has the ex- 
change constant Jrkky is expressed using the equa- 
tion [TJ We assumed that the band structure of the stud- 
ied Gei-2;-ySn3;MnyTe crystals is similar to that of GeTe. 
GeTe has got a nearly spherical valence band with max- 
imum at the L-point of the Brillouin zone^.^ the effective 
mass of the carriers in the four-fold degenerated L-band 
equal to 1.2-me [23|, where m-e is the free electron mass, 
as in the case of the GeTe crystals. The SS model as- 
sumes the finite range of the interactions by introducing 
the exponential damping factor A. In our calculations we 
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FIG. 2. The Curie- Weiss temperature normalized to the Mn 
molar fraction y as a function of the carrier concentration n 
calculated theoretically (lines) using values of the exchange in- 
tegral Jpd varying between 0.4 eV to 0.8 eV with step 0.05 eV, 
and the experimental values (points) obtained for the stud- 
ied Gei-jj-ySuajMuyTe crystals grouped with respect to the 
amount of Mn ions in the alloy y (see legend). 



assumed A to be equal to 1 nm what is a reasonable value 
for the semiconductor with relatively low carrier mobil- 
ity ^ < 100 cmV(Vs) (see Refs. US and [HI). The lattice 
parameters for the pure GeTe crystals were used in our 
calculations. 

The results of calculation of the 8(n) dependencies 
performed for values of the Mn-hole exchange integral Jpd 
varying between 0.4 and 0.8 eV are presented in Fig. [5] 
Most of the experimental points lie between the 9(n) 
curves calculated for Jpd varying between the two values 
0.4-^0.8 eV. Because of the large scatter of the exper- 
imental values we found it reasonable to focus on the 
spin-glass transition temperature Tsg estimations which 
as will be seen give more accurate results. 

The Tsg dependence on the Mn content for a few val- 
ues of the hole concentrations is presented in Fig. [31 
The experimentally determined values of Tsciv) are well 
described by the theoretical curves with the values of 
the exchange constant Jpd changing between 0.4 and 
0.5 eV. It must be noted, that for the samples for which 
the Tsg values are in the vicinity of Jpd = 0.5 eV, the 
coercive fields were the highest in the entire series of 
Gei-2:-j,Sna;Mnj,Te crystals. It may be interpreted that 
the estimated higher value of the exchange constant for 
these crystals is not an artifact but reflects their mag- 
netic properties. 

The estimated values of Jpd for the studied 
Gei-2;-j,Sna;Mnj,Te crystals are within the range be- 
tween reported in literature for Sni-2:Mn2:Te bulk crys- 
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FIG. 3. The spin glass transition temperatures Tsg as a func- 
tion of the Mn molar fraction y calculated (lines) for two dif- 
ferent values of the exchange integral Jpd and different values 
of carrier concentration n (see legend), and the experimentally 
determined values of Tsg for Gei-x-ySuxMnyTe crystals with 
different chemical content y. The points were grouped with 
respect to the different carrier concentrations (see legend). 



tals [ii] (Jpd « 0.1 eV for x<0.1) and Gei^^Mn^Te 
{Jpd Ri 0.8 ^ 0.9 eV for x > 0.1) bulk crystals [H i], and 
Jpd « 0.7 eV for x = 0.07 as obtained in Gei-2;Mn2;Te epi- 
taxial layers [2a |. 



III. 3. High field magnetization 

The magnetic properties of Gei-j;-j,Sn2;MnyTe crystals 
is the presence of high magnetic field for B < 90 kOe were 
studied at temperatures between 4.3 < T < 80 K. At mag- 
netic fields B < 600 Oe an irreversible behavior with clear 
hysteresis is observed. Typical results of hysteresis curves 
recorded for a Geo.syiSno.ogoMno.oagTe sample at a few 
temperatures are presented in Fig. [H The magnetic hys- 
teresis is not a typical feature in canonical spin-glasses 
such as CuMn [lj| but was observed for many spin-glass 
systems such as AuFe with 8% at. Fe '2&\ or NiMn with 
21% at. Mn [2^. 

The analysis of hysteresis loops showed that there ex- 
ists a correlation between the Mn content y, the values 
of spontaneous magnetization Mn, and the coercive field 
He (see Fig. HJd). When y changes between 0.012 and 
0.039 the values (measured at T w4.5 K) of Mn increase 
from 0.45 emu/g to 2.5 emu/g and He values are in- 
creasing from 80 Oe to 430 Oe. In the case of the group 
of crystals with y > 0.039 the inverse dependencies were 
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FIG. 4. a) Hysteresis loops measured at selected constant temperatures (see legend) for selected Gei-j^-ySn^^MnyTe sample 
(chemical composition shown in legend) and b) spontaneous magnetization Mr and coercive field He as a function of temper- 
ature for a few selected Gei-aj-ySnajMnyTe samples with different chemical composition (see legend). 



observed, namely, the reduction of both Mr and Hq 
with increasing y. In the case of the group of crystals 
with Mn content «0.035 and Sn content x changing be- 
tween a: = 0.090 and 0.145, the changes of the He and 
Affl can be attributed to the increase of the carrier con- 
centration. The coercive field He becomes smaller with 
increasing amount of tin ions changing from 430 Oe to 
210 Oe for crystals with a; = 0.090 and 0.145. In turn, 
the value of M^ showed the opposite trend rising from 
1.6 to 2.6 emu/g. 

The observed parameters of the magnetic hysteresis 
loop and their changes with the chemical composition of 
the alloy can be associated with the modification of the 
domain structure of the material. Furthermore, together 
with the increasing Mn content the probability for anti- 
ferromagnetic pairing of Mn ions increases. 

The He values observed in our studies are smaller 
than reported in the literature for bulk Gei-a:Mna:Te crys- 
tals [2] (iJc< 1200 Oe for a; = 0.048), but are close to 
those reported for thin Gei-a:Mna;Te epitaxial layers |25| 
(i/c<700Oefor a;<0.6). 

The typical results for the Geo.824Sno.i42Mno.o34Te 
crystal obtained in the range of magnetic fields up to 
90 kOe and temperatures lower than 80 K are presented 
in Fig. [5l The M{B) curves do not saturate and are 
nearly linear in the range of magnetic fields i? > 30 kOe. 
The lack of saturation of M reflects the fact that the 
studied alloy is not a ferromagnet but a spin glass at the 
temperatures below 50 K 

In order to model the field dependence of magneti- 
zation, we performed calculations within the molecular- 
field approximation, constructing the model described 
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FIG. 5. Isothermal magnetization curves M{B) measured 
(points) and reproduced using the mean field model described 
in the Appendix (lines) for the parameters Jpd=0.75 eV, 
y = 0.034 , p = 0.71 , A = 15 A, .lAF/ks =-20 K at different 
temperatures (see legend) for selected Gei-aj-ySn^^MuyTe sam- 
ple with a; = 0.142±0.014 and y = 0.034±0.003. 



in the Appendix. The numerical calculations were per- 
formed for the exchange constant Jpd, the Mn amount y, 
the fraction of isolated magnetic ions p, the interaction 
decay distance A, and the energy of the antiferromag- 



IV. SUMMARY AND CONCLUSIONS 
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netic interactions in the system J af specified and for 
2-10^ realizations of disorder, summing the interactions 
between the magnetic ions up to the distance of 100 A. 
The method allowed to satisfactorily reproduce the high- 
field behavior of M{H) curves, i.e. the linear increase 
of magnetization with the field. It must be emphasized 
that the method adopted is based on the molecular-field 
approximation, so that it cannot be expected to account 
for the collective excitations in the system and to yield 
the precise temperature dependence of magnetization. It 
is obvious that the low-field magnetic susceptibility is 
underestimated by the present model, which might be 
attributed to neglect of correlations. 

In Fig.|n]the isothermal magnetization curves obtained 
for the studied crystals are presented in the form of Ar- 
ret plots. With the increase of the Mn content in the 
samples the experimental M'^{B /M) curves change their 
character. Such a behavior is usually attributed to the 
change of the strength of the ferromagnetic interactions 
in the system. In the case of the Gei-a;-ySna;MnyTe sam- 
ples with Mn content y < 0.039 the observed behavior is 
similar to the one observed in many canonical spin-glass 
systems with predominance of the magnetic interactions 
with the exchange constant J > [28] . It is well known 
that in case of a spin-glass systems the M'^{B/M) curves 
on the Arrot plot lies on the bottom side of the plot what 
reflects the absence of spontaneous magnetization in the 
system. With the increase of the amount of Mn above 
2/ = 0.039 the M'^{B/M) curves on the Arrot plot shifts 
towards the top-left side of the plot. This is a signature 
that the ferromagnetic interaction dominates in the stud- 
ied samples [28]. Nevertheless, the system remains in the 
frustrated spin-glass state (see Fig.lHl). 



We presented detailed studies of the magnetic prop- 
erties of the Gei-x-ySuxM-HyTe crystals with different 
amount of both magnetic and nonmagnetic constituents 
i.e. Mn and Sn ions. 

The DC magnetometry measurements showed that the 
studied spin-glass system shows some features character- 
istic of a ferromagnetic material. Hysteresis loops with 
large coercive fields (maximum of 500 Oe) and sponta- 
neous magnetization were observed indicating that the 
studied system possessed features characteristic of a fer- 
romagnetic material. The high field magnetization curves 
showed nonsaturating behavior at magnetic fields as high 
as 90 kOe which was interpreted as a direct proof of a 
strong magnetic frustration in this system. The observed 
bifurcations of the FC/ZFC magnetization M{T) depen- 
dencies as well as the non-saturating M{B) curves are the 
features characteristic for spin-glasses. The experimen- 
tal M{T) curves were reproduced theoretically using the 
non- interacting spin wave theory. The M{B) curves were 
fitted theoretically using the molecular field approxima- 
tion. Theoretical calculations reproduced well the exper- 
imental results especially in the range of high magnetic 
fields. The linear M{B) dependence was reproduced by 
including magnetic disorder in the theoretical model and 
separate treatment of the antiferromagnetically-coupled 
pairs of magnetic ions. The proposed model gives satis- 
factory agreement with the experimental curves only in 
the case of the presence of the antiferromagnetic nearest- 
neighbor interactions in the system. This is an evidence 
that strong magnetic disorder is present in the studied 
system. 

The calculations based on a modified Sherrington- 
Southern model showed that the spin glass state is the 
preferred magnetic state for the experimentally deter- 
mined hole concentrations as well as Mn molar fractions. 
The calculations allowed us to reproduce with good ac- 
curacy the spin glass transition temperatures for all the 
studied crystals and to estimate the Mn-hole exchange 
integral value to be around Jpd = 0.45±0.05 eV for the 
studied Gei-a;-ySna;Mnj^Te mixed crystals. 
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Appendix: Theoretical model including 
antiferromagnetic interactions 

The hamiltonian of the system of interacting magnetic 
ions distributed over the lattice sites is: 



n 



\Y.j [r,, ) e. s. s, - asiiBB J2 C^ s? 



(A.l) 



^,3 



where the operators S^ describe the spins S — 5/2 dis- 
tributed over N fee lattice sites (the rhombohedral struc- 
ture of Gei-x-ySiixMuy is very close to the fee structure). 
The exchange integral is given by Eq. ([T]). The disorder in 
the site-diluted magnetic system is introduced by occupa- 
tion number operators ^i taking the values of 0, 1 for each 
lattice site i. The set of ^i values describes completely 
the given realization of the disorder in the system. 

Because of presence of the disorder, the lattice sites 
occupied by magnetic ions are inequivalent in the ther- 
modynamic meaning. The most pronounced source of 
inhomogenity is due to additional superexchange cou- 
pling occurring between the magnetic ions in nearest- 
neighbor position. In order to account for this, the sys- 
tem of magnetic impurities located on lattice sites (with 
atomic concentration equal to y) was subdivided into two 
(interacting) subsystems. For low values of y the signifi- 
cant number of impurity ions can be expected to lack the 
other magnetic impurities in the nearest-neighbor posi- 
tion. Therefore, such magnetic ions do not couple anti- 
ferromagnetically to other magnetic ions via the superex- 
change mechanism. Let us denote the total number of 
such ions (called further 'isolated magnetic ions') by Npy. 
For purely random, uncorrelated occupation of fee lattice 
sites by impurities, the probability that a selected mag- 
netic ion is isolated equals to p = {1 — y) (what yields 
the value of 0.66 for y = 0.034). However, if the distribu- 
tion of magnetic impurities is not completely random, p 
deviates from the above mentioned value {p < (1 — y) 
for ion clustering tendency, p > {1 — y) for the opposite 
tendency) . 

We assume approximately that the remaining mag- 
netic ions form isolated nearest-neighbor pairs (i.e. pairs 
of ions that do not possess further magnetic nearest 
neighbors). Within each pair, both ions are coupled via 
the superexchange mechanism, characterized by the ex- 
change integral Jaf < 0. We neglect the possibility of 

I 



formation of magnetic clusters containing more than two 
superexchange-coupled magnetic ions, which is justified 
for small y. According to this assumption, the number 
of isolated pairs is ^y{l—p). Let us emphasize that 
due to the neglect of larger magnetic clusters, the value 
of parameter p determined from the best-fit to the ex- 
perimental data may differ from the actual fraction of 
isolated magnetic ions in the real system. 

To construct the thermodynamic description of the 
model, we apply a molecular-field approximation, assum- 
ing that the total state of the system takes the form of a 
tensor product of the appropriate single-site density ma- 
trices. The quantum state of each site occupied by an 
isolated magnetic ion is described by a density matrix of 
the form 



Pj,o = exp [(Ao + gsfJ-sB) Sf //cbT] /Zq, 



(A.2) 



with Zq = J2s=-s '^■^P [('^0 + Qst^sB) s/k-oT]. The ther- 
modynamic average of the spin value of such an ion equals 
to 



(S^) = mo = SBs 



Aq + gsp.BB \ 



knT 



(A.3) 



SBs (x) being the Brillouin function for spin S. 

Let us further denote one of the (inequivalent) ions be- 
longing to an isolated nearest-neighbor pair by -f , while 
the other one by — . The corresponding density matrices 
are assumed in the form: 



p,,± = exp [(A± + gs^^BB) Sl/k^T] /Z±, 



(A.4) 



with Z± = X]s=-s'5'^P[(-^± +5sA^sB)s/fcBr]. The 
thermodynamic average of spin value of each of the ions 

is 

(Sn=-± = ^S.(^^^g^). (A.5) 

The parameters Aq, A+ and A_ are the variational pa- 
rameters of the molecular field. 

For such a model, the total entropy is a sum of con- 
tributions originating from each single isolated magnetic 
ion (— fceTr [pi,o In p^.o]) and from each isolated magnetic 
pair (-fceTr [pj_+lnpj_+] - fceTr [pj__ lnpi^_]). The to- 
tal entropy of the system is: 



S — k^ 2_^ ii 



l-p 
plnZifi H -— (lnZj,+ + lnZ,;._) + 



Ao + gsfJ'BB 
P k^T "^" 



1 -p fX + + gsfJ'BB 
ksT ' 

I 



A_ + gsfiBB 
knT 



(A.' 



The enthalpy of the model is a thermodynamic average interaction between the isolated ions, between the iso- 
of the Hamiltonian, containing the terms describing the lated magnetic pairs as well as the interaction between 



isolated magnetic ions and isolated pairs. The total en- thalpy amounts to: 
I 



("H) = -^Cj < ■^P^ml^^jJ{R,j) +pmogsfJ-BB + -(1 -p)JAFm+m- 

i K j 

H ^gsfJ'BB{m,+ + m^) + p——mo{m+ + m_) ^ ^^ 7(7?,^ ) + 

i 



(A.7) 



Let us note that in the formula above, as well as in the 
rest of the Appendix, the nearest-neighbors of the lattice 
site i will be excluded from summation over j. 

The total Gibbs free energy of the system is obtained 
from the expression G — (H) — ST. The variational min- 
imization of the Gibbs energy with respect to molecular 
field parameters yields the set of coupled self-consistent 
equation of the form: 



j 
A_ = JAFm+ + pmo ^ ^j J{Rj), 



(A.9) 



Ee^ Ao 



pmo 



1-p 



(to_| 



which need to be solved together with the conditions 
(|A.3p and (IA.5I) . The solution allows us to calculate the 
total magnetization in the system for a given realization 
^^jJ{Rij) /• = of disorder as: 



E ^M A+ - JApm^ - pniQ ^ £,jJ{Ri 



M ■ 



nmgSfJ-B 



pmo 



{m+ + m_) 



(A.IO) 



E ^M A- - JAFm+ - pmo ^ Cj J [Rij ) > = 0. 



Then we are looking for the solution when all the terms 
under the sum over i vanish simultaneously, i.e. we ob- 
tain the following set of equations : 



Ao 



pmo 



1-p 



(m+ + TO_) 



^JO-^l^jO 



(A. 8) Finally, to obtain the disorder-averaged site mag- 
netization, we average the values obtained from the 
Eq. ()A.10|) for a sufficiently large number of realizations 
of disorder. Each realization of disorder is simulated 
numerically, by allowing the parameters ^j to take the 
values of 0, 1 according to the probability distribution 
P i^j) = yS iS,j ~ 1) + (1 ^ Z/) (^ (^j): with the total con- 
centration of magnetic impurities equal to y. 
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